arXiv:1504.05356v2 [gr-qc] 12Jun2015 


Thermal corpuscular black holes 

Roberto Casadio“^* Andrea Giugno“^| and Alessio Orlandi“^^ 

Dipartimento di Fisica e Astronomia, Universitd di Bologna 
via Irnerio 46, 1-40126 Bologna, Italy 

^I.N.F.N., Sezione di Bologna, 
via B. Pichat 6/2, 1-40127 Bologna, Italy 

June 15, 2015 


Abstract 

We study the corpuscular model of an evaporating black hole consisting of a specific quantum 
state for a large number N of self-confined bosons. The single-particle spectrum contains a 
discrete ground state of energy m (corresponding to toy gravitons forming the black hole), and 
a gapless continuous spectrum (to accommodate for the Hawking radiation with energy uj > m). 
Each constituent is in a superposition of the ground state and a Planckian distribution at the 
expected Hawking temperature in the continuum. We first find that, assuming the Hawking 
radiation is the leading effect of the internal scatterings, the corresponding A-particle state 
can be collectively described by a single-particle wave-function given by a superposition of a 
total ground state with energy M = N m and a Planckian distribution for E > M at the same 
Hawking temperature. From this collective state, we compute the partition function and obtain 
an entropy which reproduces the usual area law with a logarithmic correction precisely related 
with the Hawking component. By means of the horizon wave-function for the system, we finally 
show the backreaction of modes with uj > m reduces the Hawking flux. Both corrections, to the 
entropy and to the Hawking flux, suggest the evaporation properly stops for vanishing mass, if 
the black hole is in this particular quantum state. 


1 Introduction 

Recently Dvali and Gomez [1] proposed the idea that a black hole can be modelled as a Bose- 
Einstein condensate (BEG) of marginally bound, self-interacting gravitons. In fact, these bosons 
can superpose in a single small region of space, effectively giving rise to a gravitational well, whose 
depth is proportional to the total number of constituents. Furthermore, in the mean field picture, the 
Bogoliubov modes that become degenerate and nearly gapless, represent the holographic quantum 
degrees of freedom responsible for the black hole entropy [2] and the information storage [3]. Among 
the merits of this model, there is thus the fact that it provides the new perspective of describing the 
Hawking radiation as an emission of quanta which are already present in the system, rather than 
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being created out of the semiclassical vacuum. The result in Refs. [1, 3] have thus triggered a number 
of developments [4, 5, 6, 7, 8, 9] and possible cosmological implications were also investigated in 
Refs. [10, 11, 12, 13]. However, the lack of a geometrical description of the black hole space-time 
makes it difficult to assess the presence of a horizon in this model. 

A possible way to make contact with the usual point-of-view of general relativity was later 
proposed [8], using the formalism of Refs. [14]. This approach can in principle be applied to any 
quantum state and introduces the idea of a “fuzzy” gravitational radius i?H the quantum level, 
starting from the Einstein equation that determines the classical Misner-Sharp mass for spherically 
symmetric systems and therefore determines the location of trapping surfaces. For example, if a 
particle is in superposition of many energy eigenstates, to each eigenstate there will correspond a 
different value of i?H) with a probability amplitude given by the corresponding spectral coefficient. 
The particle’s size could then be smaller than the mean value of i?H with finite probability. Simple 
evaluations showed that only particles with a mass of the order of the Planck scale or larger, are 
likely to appear in a black-hole state. Moreover, for astrophysical masses, the relative uncertainty 
in the horizon size can be acceptably small for semiclassical black holes provided the system is made 
of a large number N of very light condensed bosons [8]. In fact, in this case the horizon relative 
uncertainty decreases with N, which supports the idea that large BECs of gravitons at the critical 
point can be viewed as semiclassical black holes in the large-limit. 

More specifically, in Ref. [8], we considered two possible candidates for the quantum state 
representing a black hole and its Hawking radiation: one in which the Hawking quanta have a 
Gaussian distribution and one with Boltzmann distribution in energy. Unlike the former, the latter 
state explicitly contains the Hawking temperature, but cannot be used beyond perturbation theory 
(see Appendix A). Hence, in this work we shall consider yet a different quantum state for the 
BEG and Hawking radiation, namely we shall assume each component in the BEG has a Planckian 
probability amplitude to be in an excited state. This choice does not suffer from the drawback of 
the Boltzmann distribution, and will allow us to connect nicely with the standard thermodynamic 
picture, and predict the vanishing of the specific heat for small (possibly vanishing) black hole 
mass. Finally, it will also let us estimate the backreaction of the emitted quanta by determining 
the horizon wave-function and mean gravitational radius of the system, and show that the Hawking 
flux vanishes at (vanishingly) small mass, as one should expect from energy conservation. 

We shall start by briefly reviewing the original model of Refs. [1] in the next section, and then 
build a specific quantum Ai-particle state that could refine their results in section 3. In fact, we 
shall show that such a state entails the Bekenstein-Hawking entropy with a logarithmic correction 
and Hawking flux of outgoing particles in section 4. Section 5 is instead devoted to the calculation 
of the expectation value of the gravitational radius of the system, and the (small) deviation from 
the expected classical values will be used to estimate the backreaction of the Hawking flux. Finally, 
we will conclude with some comments in section 6. 

2 Black holes as BECs 

Before we describe the specific quantum state we shall use to model an evaporating black hole, let 
us review the basics of the model of Refs. [1]. The whole construction starts from the assumption 
that the gravitational interaction, approximated by the Newtonian potential 

UN(r) , (2.1) 

r 
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should be strong enough to conhne the gravitons themselves inside a hnite volume. If so, the 
gravitons will acquire an effective mass m related to their characteristic quantum-mechanical size via 
the Compton/de Broglie wavelength Am — him = (.p mp/m. In fact, N gravitons could superpose in 
a spherical volume of approximate radius Am, and total energy M = N m. Within this Newtonian 
approximation, the effective gravitational coupling constant is thus given by a = |I4^(Am)|/.N ~ 
= m^/mp, and the average potential energy per graviton can be estimated as 

U ~ mVN(Am) — —Nam . (2.2) 

If we hnally assume the gravitons are “marginally bound”, so that 

Ek ^ 0 , (2-3) 

where Ek — m is the graviton kinetic energy, we obtain the “maximal packing” condition 


N" a ~ 1 , 

which imply that the effective graviton mass and total mass of the black hole scale like 

m ~ —JL and M = N m ^ \/]V m^ ■ 

Vn ^ 

The horizon’s size, namely the Schwarzschild radius 

Rr = 2 ip - , 

mn 


(2.4) 

(2.5) 

( 2 . 6 ) 


is therefore of the order of the de Broglie length Am — ip mpjm of the gravitons and is quantised 
as commonly expected [2], that is 

Rn-VNip. (2.7) 

For our purpose, it is of particular importance that the Hawking radiation and its negative 
specihc heat spontaneously result from the quantum depletion of the condensate for the states 
satisfying Eq. (2.3). At first order, reciprocal 2—7-2 scatterings inside the condensate will give rise 
to a depletion rate 


1 




1 


( 2 . 8 ) 


VNip ’ 

where the factor N~‘^ ~ the second factor is combinatoric (there are about N gravitons scattering 
with other — 1 ~ gravitons), and the last factor comes from the characteristic energy of the 
process AE ~ m. The amount of gravitons in the condensate will then decrease according to [1] 

iV~-r~- j^ + 0(N-^). (2.9) 

VNip 

As explained in Refs. [1], this emission of gravitons reproduces the purely gravitational part of the 
Hawking radiation and contributes to the shrinking of the black hole according to the standard 
results 

.3 


M ~ mr 


N 


■\/N j v ^p 


rup 

' NL 




£pM2 ■ 


From this flux one can read off the “effective” Hawking temperature 

,2 


Th ~ 


mr 


mr 


~ m ~ 


SttM y/N ’ 

where the last expression is precisely the approximate value we shall use throughout. 


( 2 . 10 ) 


( 2 . 11 ) 
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3 BEC with thermal quantum hair 

Like in Ref. [ 8 ], we start by considering a system of N scalar particles, i = 1,... ,N, whose dynamics 
is determined by a Hamiltonian Hi. We do not need to specify the latter, and will refer to such bosons 
as the “toy gravitons”. Since we want our particles to be marginally bound according to Eq. (2.3), 
we assume the single-particle Hilbert space contains the discrete ground state \m), defined by 

Hi\m) = m\m) , (3-1) 

and a gapless continuous spectrum of energy eigenstates |wj), such that 

\^i) — 1^^) ? (3'2) 


with uji > m. This continuous spectrum is meant to reproduce the depleted toy gravitons that will 
escape the BEC. Each particle is then assumed to be in a state given by a superposition of \m) and 
the continuous spectrum, namely 

,r,(i)> _ l™> +71 

where 71 > 0 is a real parameter that weights the relative probability amplitude for each particle 
to be in the continuum rather than ground state. 

In Ref. [ 8 ], we considered both a Gaussian and a Boltzmann distribution for the continuum. Since 
the Boltzmann distribution leads to inconsistent results for general values of 71 (see Appendix A), 
we instead assume a Planckian distribution at the temperature Th = m [from Eq. (2.11)] for the 
continuum part, namely 


(3.3) 




Nh 

™3/2 


do;,- 


{oji - m) 


roo 

Hh d£iG(Li) |£i) , 
Jo 




m 

'OQ 


(3.4) 


where the normalisation ATh = — 6 C(3) — 1.06, we introduced the dimensionless variables 


£,: = 


w,- — m 


m 


and the states |£j) = Iwj), such that the identity in the continuum can be written as 


POO POO 

/ d£i|£i)(£i|= / dcjj |wi) (cjil = I 

^0 J m 


(3.5) 


(3.6) 


Since we assume that this “thermal hair” arises because of the scatterings among the scalars inside 
the BEC, the parameter 71 should be related to the toy gravitons self-coupling a, and therefore 
vanish for a —)• 0 or A^ —)• 0. This in turn would mean that in principle 71 should also depend on 
N, for example if the maximal packing condition (2.4) remains satisfied along the evaporation (see 
Ref. [9] for a more detailed study of the scatterings inside the BEC). We however prefer to keep 
it as a free variable in this work, so that one could relate it a posteriori to known features of the 
Hawking radiation, at least in the large N limit. 
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In the following, we shall just assume that 71 7 ^ 0 yields a sufficiently good approximation of 
the leading effects due to these bosons self-interactions and therefore treat the BEC as made of 
otherwise free scalars. The total wave-function of the system of N such bosons will correspondingly 
be approximated by the totally symmetrised product 


I^Af) ^ 


1 

M 


N 


N 

E 

Li=i 


1 ^ 




(3.7) 


where Ylfa} denotes the sum over all of the A^! permutations {cTj} of the N terms inside the square 
brackets. Upon expanding in powers of 71 , we obtain the exact expression 


N 


1'^^)"" mE 

■ ffid Li=l 
N 


N 


m) 




N 


77 


{a,} U=2 

N 


ffid Li=3 


poo 

\m)(E) / d 8 iG( 8 i)|£i) 

Jo 

poo poo 

\m)(E) / d£iG(£i)|£i)® / d £2 ^(£ 2 ) 

Jo Jo 


+ ... 






^ poo 

(g)|m)(g)/ d£,G(£,)|£, 

=j+i 3=1 


+ ... 

+7f ^ 

fV! 


N V N 


V /*oo 

fE ® / d£,G(£,)|£,) 

• {<xa Li=i -^0 


(3.8) 


where we omitted the overall normalisation constant of 1/(1 -|- 7 i)'^^^ for the sake of simplicity. 

Since we are effectively including the interaction into terms proportional to powers of 71 , the 
spectral decomposition of this A^-particle state can be obtained by defining the total Hamiltonian 
simply as the sum of N single-particle Hamiltonians, 


N 

H = ^H,. (3.9) 

i=l 

The corresponding eigenvector for the discrete ground state with M = N m is thus defined by 


H\M) = M\M) , (3.10) 

and the eigenvectors for the continuum by 

H\E)=E\E). (3.11) 


The spectral coefficients are then computed by projecting ITtv) on these eigenvectors. For E = 
M = N m, again neglecting an overall normalisation factor, we obtain 


N 


C(M) c, T (Ml 

Wi} 


N 


\m) 


Li=l 


= 1 


(3.12) 
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as expected. For energies above the ground state, E > M = N m, they are likewise given by 

C{E > M) = {E \^n) , (3.13) 

and it is again convenient to employ the dimensionless variables (3.5), and further introduce 

E-M 


£ = 


m 


(3.14) 


which lead to 

C(£ > 0) ~ 7iXhG(£) 

roo 

WiK G(£i)G(£-£i)d£i 

Jo 


+ ... 


roo roo / N \ 

Gi8,i)d8,ix---x G(£^)d£^<5 £-V£i 

Jo Jo \ , 

N 

E^rc«(£). 


(3.15) 


n=l 


where all the coefficients in this expansion can be written as 

roo roo / \ 

Cn = y^J^ G(£i)d£i X ••• X G(£„_i)d£„_iGf£-^£j . 

We then note that each integral in £j is peaked around £j = 0, so that we can approximate 


(3.16) 


n—1 


£ - £ 


a £-E£. 


\ i=i / exp 
for 2 < n < A'". We thus find 


exp{£}-l 


= G{E) , 


(3.17) 


TT 


2 \ n-l 


Upon rescaling 


~ Nh ( Nh y ) G(£) = (1.75)--1 Xh G(£) . 


N 

7-0.57 ^( 1 . 7571 )^' , 
i=i 


(3.18) 


(3.19) 


and switching back to dimensionful variables, we immediately obtain 

Nh {E-M)/m 


G{E > M) - 7 


m exp {{E — M)/m} — 1 


(3.20) 


This approximation was checked numerically to work extremely well for a wide range of N (see 
Appendix B). 
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The result (3.20) means that we can collectively describe the quantum state of our A^-particle 
system as the (normalised) single-particle state 


\M) +'f\'ip) 


(3.21) 


where 


^ /■- {E-M)/m 

exp{(E-M)/m}-l' ^ 


(3.22) 


In other words, from the energetic point of view, the BEC black hole effectively looks like one 
particle of very large mass M = N m in a superposition of states with “Planckian hair”. If we 
wanted to have most of the toy gravitons in the ground state, we could just assume that 7 <C 1 . 
However, nothing prevents one from also considering the above approximate wave-function for 7 ~ 1 
or even larger. 


4 Partition function and entropy 

It is now straightforward to employ the effective wave-function (3.21) to compute expectations 
values, such as 

. -T_/M+ 7 ^r[ i^d 4 

l + 7^\ m Jm [exp {(.E - M)/m} - ij j 

if / pcxD roo \ 

= M + G^{e,)de. + mJ 0^(8) £d8j 


1 + 7' 


^ M(1+ 7 ^)+ ( 6C(3)- 




where we used ( 1 + 7 ^)“^ ~ and discarded higher powers of the parameter 7 (see Appendix C 

for the integrals). 

Since 

1 + 72 I m Jm [exp {(E; - M)/m} - ij j 

-I r / roo poo poo \ 

= + G'2(£)d£+ 2Mmy {E) E dS. + J G^{E) dS, j 

= mlN 1 + 2^1^^ (6C(3)-^) (7r"-90C(5)) + 0 ( 7 "), (4.2) 

it is easy to obtain the uncertainty 

AE = ^/{HY-{Hf = 7 ^ :N-h (vr^ - 90C(5)) - (6C(3) - + 0 ( 7 ^) 

7T) 

~ 0 . 767 ^. (4.3) 

Wn 


7 



(4.4) 


This allows us to recover the large N result [ 8 ] 


AE 

W) 


~ ^ + 0(7^) . 


Adopting a thermodynamical point of view, we can use (4.1) to estimate the partition function 
of the system according to 


(H) 


d 


logZ(/3) , 


(4.5) 


where (3 = ^ = m ^ ~ ^/N/nip. We then have 


{HW)) 






d 


^771^/3^+A7Mog(/c/3) , 


(4.6) 


where A = [6(((3) — 7r^/15] ~ 0.81, and k is an integration constant with the dimensions of a 

mass. It is then straightforward to obtain 


log Z(/3) = -^m^/3^ -Aj^ log(/3 k) , 

(4.7) 

and, if we simply set k = nip, we see that 


Z{(3) = (mp/3)-^^' , 

(4.8) 


contains two factors. One goes like 1//3 = Th at some power, like the thermal wave-length of a 
particle does, and it is given by the contribution of the excited modes to the energy spectrum. This 
is consistent with the fact that such modes propagate freely, since we did not associate an external 
potential to our model ^. The exponential damping factor is instead due to the presence of bound 
states localised within the classical Schwarzschild radius i?H- 

It is now possible to compute the statistical canonical entropy 


5 (/ 3 ) = /32 


dF{/i) 

d/3 


(4.9) 


where F(/3) is the Helmoltz free energy F = —(1//3) logZ. It is straightforward to get 

S(/3) = ^m^/3^-Aj^log(mp/3)+Aj^ , (4.10) 

which is nothing but the usual Bekenstein-Hawking formula [2] plus a correction scaling like a 
logarithm. In fact, we can write /? ~ m~^ as a function of the Schwarzschild radius. 


/? 


2£pmp 


rescale the constant and define the area of the horizon as Ah = dvri?^, hence yielding 


(4.11) 


S 



A 


7 ^ log 


Ah \ 
16 y 


1 


In other words, we are approximating the grey-body factors for bosonic Hawking quanta with one. 


(4.12) 



Let us conclude with a few remarks. First off, the final expression of the entropy depends 
on the “collective” parameter 7 defined in Eq. (3.19), rather than the relative probability 71 for 
each constituent to be a Hawking mode, the latter being in turn determined by the details of the 
scatterings that occur inside the BEC. One could thus speculate that, even if we were able to 
detect the Hawking radiation, the details of these interactions would hardly show directly. The fact 
remains, though, that 71 = 0 implies that also 7 = 0 , and and the logarithmic correction would 
therefore switch off if the scatterings inside the BEC were negligible (and the Hawking radiation 
were absent). Finally, we stress again that, in this corpuscular model, the total number N of bosons 
is conserved, since both the black hole and its Hawking radiation are made of the same kind of 
particles. In connection with this, we notice that the specific heat is in our case given by 


Cv 


d{H) 

dT 




d{H] 

df3 


~ —m. 


, + yi 7^ 


(4.13) 


which is negative for small 7 and large [3 (or, equivalently, large N), in agreement with the usual 
properties of the Hawking radiation, but vanishes for j5 ~ that is for Nc ~ 7 ^. If, for 

simplicity, we assume 1.75 71 ~ 1, so that 7 ~ iV from Eq. (3.19), we thus obtain the specific heat 
vanishes for iVc ~ 1, as one would naively expect the Hawking radiation switches off when there are 
no more particles to emit. In fact, this result qualitatively agrees with the microcanonical picture 
of the Hawking evaporation ^ and with the estimate of the backreaction in the next section. 


5 Horizon wave-function and backreaction 

We can now investigate the presence of a trapping surface in our quantum state by means of the 
horizon wave-function formalism. Since the present case is very similar to the ones considered in 
Ref. [ 8 ], we refer the reader to that paper for more details (see also Refs. [15] for a similar picture). 

The main idea behind that formalism, is to consider the relation (2.6) that defines the gravita¬ 
tional radius of a spherically symmetric system as an operator equation, that is 

fn = 2£p — . (5.1) 

rrip 

Note that the classical gravitational radius rn = rn^r) = 2ipM(r)/mp can be introduced for 
any spherically symmetric system with mass function M = M{r) [16], and it does not necessarily 
represent the size of a horizon. In fact, it only represents the possible location of a trapping surface 
if r}i{r) = r for some value of the areal coordinate r. The above definition allows one to introduce 
a wave-function for the gravitational radius simply given by 

TH(rH) ^ C{mprn/2£p) , (5.2) 

where C is precisely the spectral coefficient we computed in the previous section, and the normali¬ 
sation will be fixed assuming the Schrodinger scalar product 

roo 

(Th I 4>h) = dvr / TH(rH)d»H(rH)?’HdrH • (5.3) 

Jo 

^For the details, see Refs. [17] and references therein, where the black hole and its Hawking radiation where also 
assumed to be of the same nature, to wit strings. 
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It is important to recall that the main result in Ref. [8] was that the gravitational radius associated 
with the quantum iV-particle states considered therein indeed has a size very close to the expected 
classical value R}i, with quantum fluctuations that die out for increasing N, namely 

(rn) + , (5.4) 

as one should expect in the context of semiclassical gravity. Since the source represented by such 
quantum states is by construction (mostly) localised within R}i, this shows that r = (fn) is a 
horizon for the system. 

Given the spectral coefficients (3.20), the corresponding horizon wave-function reads 

|4'h) = - , „ , (5.5) 

y'l+T 

where the state |i?H) represents the discrete part of the spectrum and is defined so that 

{Rnl rn |i?H) = Rk , (5.6) 


and the states 

ipnirn > Ru) = (m | V'h ) = 


N {m - RR)/{ 2 £pm/mp) 


4 TT R^ gjjp 


f 1 

2£p m/rrip J 


- 1 


(5.7) 


which describe the contribution from the excited Hawking modes. As usual 'k(rH < Rn) — 0 and 
the normalization is fixed by the scalar product (5.3). The expectation value of the gravitational 
radius is now easily calculated, 


roo 

(fn) = dvr / |4'H(rH)|^ 

JRh 


1 


Rh + 7^ ^ 


1+7^ 


P 


= Rh 

= 


1 + 


7 ^ 3X^ 


1-^72 N 


P 


eP -1 
2 


{p + Nfdp 


37 ' 


eP -I 
4 


pdp 




i + ^:nM6C(3)--)|+o(7^), 


TT 

I5 


where we defined 


P = 


rn - -Rh 


2 ip m/rup 

together with the relation (2.6). We thus see that 

{fn) - Rn^^l^ ^ + O 


N 


In the same way we can compute 


{fl) = R? 


H 


TT 


l + 4yNA(6C(3)-- - 


1 




(5.8) 


(5.9) 


(5.10) 


(5.11) 
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which leads to 


Arn 





7^H 

Vn 



TT^ 


6C(3)-- 


+ 0 



(5.12) 


and, omitting a factor of order one, 


ArH 

(^h) 




1_ 

N 


+ 0 ( 7 '), 


(5.13) 


like we obtained in Ref. [ 8 ]. 

We can now note that (rn ) > albeit bny a small amonnt for large A, which is in agreement 
with the fact that Hawking qnanta mnst also contribnte to the total gravitational radins. The 
constitnents will therefore be bonnd within a sphere of radins slightly larger than the pnre BEC 
valne Rhi and the typical energy of the reciprocal 2 — 7-2 scatterings will correspondingly be smaller. 
This gives rise to a slightly smaller depletion rate, namely 


r ~ 


1 ^ 1 


1-^^h(6C(3) 



(5.14) 


Note that, if we were to trnst the above relation for small A, the flnx wonld stop for a nnmber of 
qnanta 


Nr. 


37 ' 


6C(3) 



~ 2.47^ . 


(5.15) 


If we farther recall the relation (3.19), and again assnme for simplicity 1.75 71 < 1, so that 7 < IV, 
we obtain W > 1 , which leaves open the possibility that the flnx stops for a finite nnmber of 
constitnents. Of conrse, it is not gnaranteed the above approximations still hold for snch small 
valnes of A, bnt we recall that it is a general resnlt of the microcanonical description of the Hawking 
radiation that the emitted flnx vanishes for vanishingly small black hole mass [17]. Fnrthermore, 
this behavionr wonld be in agreement with the thermodynamical analysis of the previons section, 
and the vanishing of the specific heat for A = W ~ 1. In any case, Nc depends on the collective 
parameter 7 , althongh the latter is more directly related with the single-particle 71 for small N. 


6 Conclusions and outlook 

We have analysed a candidate qnantnm state to represent an evaporating black hole made of N toy 
scalar gravitons, along the lines of the BEC pictnre pnt forward in Refs. [1], and improving on the 
cases previonsly considered in Ref. [ 8 ]. In particnlar, we fonnd that a Planckian distribntion at the 
Hawking temperatnre for the depleted modes allows one to recover the main known featnres of the 
Hawking radiation, like the Bekenstein-Hawking entropy, the negative specific heat and Hawking 
flnx, for large black hole mass, or, eqnivalently, large N. 

Moreover, one also obtains that the above A-particle state can be collectively described very 
reliably by a one-particle wave-fnnction in energy space. From this approximate description, we 
obtained leading order corrections to the energy of the system that give rise to a logarithmic con- 
tribntion to the entropy. This contribntion wonld ensnre a vanishing specific heat for N of order 
one, when we expect the evaporation stops. In fact, this qnalitative resnlt was farther snpported by 
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estimating the backreaction of the Hawking radiation onto the horizon size by means of the horizon 
wave-function of the system. Since the Hawking flux depends on the energy scale of the BEC [1], 
which is in turn simply connected to the horizon radius, the extra contribution to the latter by 
the depleted scalars precisely reduces the emission. Upon extrapolating the picture down to small 
values of N, this reduction should eventually stop the Hawking radiation completely. We think it 
is important to emphasise that this is exactly what one expects from the conservation of the total 
energy of the black hole system [17]. 

We would like to conclude by remarking that the analysis presented here does not explicitly 
consider the time evolution of the system or, in other words, it provides a possible picture of the 
black hole and its Hawking radiation at a given instant in time. One could however conjecture that, 
as long as the Hawking flux is relatively small, a reliable approximation of the time evolution will 
be simply given by the standard equation (2.10) for large mass and N, or by the correspondingly 
modified expression that follows from Eq. (5.14) for N approaching one. Our analysis however does 
not make use of any specific determination of the coefficients 71 and 7 in terms of the details of 
the scatterings that give rise to the Hawking radiation (and we therefore assumed no dependence of 
these parameters from the number N). A very interesting attempt at understanding the Hawking 
radiation directly from the two- and three-body decays inside a BEC recently appeared [9], where 
however no general relativistic effect is included. The estimate for the flux we provide in section 5 
could therefore be viewed as complementary to that in Ref. [9], and both approaches in fact agree in 
predicting a vanishing decay rate for vanishing number of constituents. Moreover, we suspect that, 
for very small N, the BEC model could reduce to the description of black holes as single particle 
states investigated in Refs. [4, 14] and [18]. A detailed analysis of the possible transition between 
the two regimes is left for future investigations. 
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A Boltzmann distribution 


It is easy to see that the Boltzmann distribution considered in Ref. [8] can only be employed in the 
perturbative limit 71 <C 1. In fact, if one replaces the Planckian distribution 

G(8 )^e-^ (A.l) 


in Eq. (3.4), and expands the corresponding C{E > M) in powers of 71 like in Eq. (3.15), the 
analogous coefficients Cn will diverge when n > 2, namely 


COO coo coo / ^ \ 

Cn o: J d8ie-^i J d 8 . 2 e-^^ x... x J da^e"^" 5 f £ - ^ 8 * J 

COO 

= e~^ / d£i —)■ 00 . 


(A.2) 


i=i -^0 


It then follows that, for 71 ~ 1, and the probability for each constituent to be a Hawking mode is 
not small, the Boltzmann approximation should indeed be replaced by the Planckian distribution 
employed in this work. 
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(a) Cs 


(b) Cio 


Figure 1: Numerical estimate of some coefficients in the expansion (3.15) (black lines), and their 
analytical approximation (red lines), for = 10 (energies are in Planck units). 

B Monte Carlo estimate of spectral coefficients 

We have checked the analytical approximation (3.20) by computing directly some of the spectral 
coefficients (3.15) by means of a Monte Carlo algorithm. Fig. 1 shows the numerical estimates of 
C 5 = C^{E — M) and Cio = Ciq{E — M) along with their analytical approximation for N = 10. 
This comparison immediately shows that the coefficients Cn indeed have a “Planckian” shape, in 
agreement with their analytical approximation. Finally, Fig. 2 shows the whole coefficient C{E > 
M) for N = 10. The same quantities for larger values of N, up to = 100 have also been computed 


C(E>M) 



Figure 2: Numerical estimate of C{E > M) (black line), and its analytical approximation (3.20) 
(red line) for A^ = 10 (energies are in Planck units). 

and the results are qualitatively the same. 
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C Useful integrals 


In this work we made use of the following results: 


JO 

roo 


X 


/o 

(■oa 


e* - 1 

X 


- 1 


- 1 


1 dx = 

I[7r2-6C(3)] , 

(G.l) 

2 

xdx = 

6C(3)-^. 

(G.2) 

x^ dx = 

4p‘-90C(5)] . 

(G.3) 
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